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A dual formulation of group field theories, obtained by a Fourier transform mapping functions 
on a group to functions on its Lie algebra, has been proposed recently. In the case of the Ooguri 
model for SO(4) BF theory, the variables of the dual field variables are thus so (4) bivectors, which 
have a direct interpretation as the discrete B variables. Here we study a modification of the model 
by means of a constraint operator implementing the simplicity of the bivectors, in such a way that 
projected fields describe metric tetrahedra. This involves a extension of the usual GFT framework, 
where boundary operators are labelled by projected spin network states. By construction, the 
Feynman amplitudes are simplicial path integrals for constrained BF theory. We show that the 
spin foam formulation of these amplitudes corresponds to a variant of the Barrett-Crane model for 
quantum gravity. We then re-examin the arguments against the Barrett-Crane model(s), in light of 
our construction. 



Introduction 

Group field theories (GFT) represent a second quantized framework for both spin networks and simplicial 
geometry |fj, |7|, |9( , field theories on group manifolds (or Lie algebras) producing Feynman amplitudes which 
can be equivalently expressed as simplicial gravity path integrals @ or spin foam models [ioj ]. in turn 
covariant formulations of spin networks dynamics [l| 3 . Most spin foam (and GFT) models for 4d gravity 
are obtained starting from models describing topological BF theory and adding constraints on the 2-form 
B field, imposing it comes from a tetrad 1-form as B = *e A e. This is based on the Plebanski formulation 
[H O of General Relativity as a constrained BF theory. A classically equivalent modification of the same 
action is obtained by adding a topological term to it, which vanishes on-shell, with a coupling constant 
named Immirzi parameter [111 ] . This gives the so-called Hoist action, which is the classical starting point 
of loop quantum gravity (LQG)[1(. The construction of spin foam amplitudes takes place in a simplicial 
context, by first assigning to the 2-dimensional faces of a simplicial complex data representing the discrete 
analogue of the fields of BF theory, that is Lie algebra elements representing the B field, and group elements 
defining a discrete connection, and then defining a quantum amplitude imposing the Plebanski constraints 
on them. Ideally, one would want to constrain the discrete data at the classical level and then quantize the 
resulting geometric structures, as one would do in a simplicial path integral. However, this has proven very 
difficult up to now, and the usual route [3, Qjl H3 starts with a quantization of topological BF theory and 
proceeds by imposing the constraints at the quantum level, that is at the level of quantum states, hoping to 
obtain the correct result nevertheless, and despite the ambiguities involved in this way of proceeding. More 
precisely, the existing models [22| (which exist in both euclidean and minkowskian signatures) result from 
the imposition of constraints via the solution of operator equations (e.g. via a master constraint technique) 
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3 They can also be seen, from a more statistical field theory perspective, as a higher dimensional generalization of matrix 
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or from constraining the quantum labels of BF coherent states, identified with the quantum analogue of 
the bivector fields of this theory, and are further distinguished by the presence and value of the Immirzi 
parameter. In absence of the Immirzi parameter 7 (on in the 7 = 00 sector of Hoist theory), one gets the 



Barrett-Crane model [17[ by using the operator method, and the so-called FK model 1221 via the coherent 
state method. With finite 7 > 1 the operator method gives the so-called EPRL model 22] or an extension of 
the FK model. These two models remarkably coincide for 7 < 1, and in particular for 7 = (corresponding 
to the 'topological sector'of Hoist gravity). The new models thus succeed in incorporating the Immirzi 
arameter in the spin foam (and GFT pjj [2t| HH) formalism, attempts at which had started very early 
37, 38]. 



In parallel with the above results, the simplicial path integral route has also been explored, providing 
important insights [2l|, |3(3|. In particular, the analysis of [2l[ clarified to a great extent the (quantum) 
simplicial geometry of the Barrett-Crane model. The results we present in this paper can be seen as a new 
take on and a generalization of the results presented there. A decisive impulse to this line of research has 
been given recently by the definition of a metric representation of GFTs H 1 , in terms of B variables, and the 
proof of the exact duality of simplicial path integrals and spin foam models, within a GFT context, based 
on non-commutative tools [13Ml6| . then applied also to LQG states [31] . Despite being phrased in terms of 
simplicial geometry and in seemingly classical terms, as appropriate for a path integral quantization, also 
this method of quantization necessarily relies on a certain choice of quantization map (which dictates, for 
example, choices of operator orderings. The one that seems [l5j to be at the root of the non-commutative 
Fourier transform and star product on which these results rely is the Duflo quantization map jS^. This 
map is the most mathematically natural one for systems whose classical phase space is (based on) the 
cotangent bundle of some group manifold [36j j , has been successfully applied to the path integral formulation 
of Chern-Simons theory |35| . and argued to be relevant also for loop quantum gravity |3q ]. Further support 
for the interpretation and use of these tools comes from their application to the much simpler case of the 
quantization of a particle on the sphere, where they allow to derive a complete and correct canonical path 
integral formulation of the dynamics [331] 2 . 

The EPRL and FK models [22| have now replaced the BC model in the interest of the spin foam community. 
On the one hand, this is due to the mentioned success in incorporating the Immirzi parameter in the 
formalism, and to the consequent close relation between the boundary states associated to these new models 
and those of canonical LQG. On the other hand, this is due to a variety of results and arguments that 
have been put forward as implying the failure of the Barrett-Crane model in correctly describing quantum 
geometry, and thus in representing a compelling candidate for a quantum gravity model 3 . 
These (numerous) arguments (related to each other) concern: the boundary states of the model (used also 
to define the vertex amplitude) and the quantization procedure used to obtain such states 22], the role 
of degenerate geometries [26| . the possible lack of (discrete analogues of) needed secondary second class 
constraints [39|, |40j, and, most importantly, the encoding of simplicial geometry in the quantum amplitudes. 

This last issue, often referred to as the 'ultralocality of the BC model' can be articulated as: 

a) 4-simplices speak only through face representations, i.e. triangle areas |25j |: 

b) bivectors associated to same triangle in different simplices are not identified [U [22J ; 

c) normal vectors to the same tetahedron seen in different 4-simplices are uncorrelated 0, |43| ; 

d) the simplicity constraints are imposed in a non-covariant fashion 

e) because of this non-covariance, there are missing constraints over the connection variables [f| . 



1 For earlier attempts aiming at the same result, see |24l . 

2 This metric formulation also allows to identify diffeomorphism transformations at the GFT level, and to link nicely their 
features across the canonical, spin foam and simplicial gravity formulation of the same theory [3^1 . 

3 These arguments, mostly, refer to the BC vertex in the spin representation, and to the choice of data entering it, also because 
there is no definite consensus on the other ingredients entering the model, i.e. 'measure' factors. 
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The last two issues had been identified in [5| , in a first application of the non-commutative metric formalism 
to GFT gravity models, which has the advantage of making the simplicial geometry behind spin foam models 
manifest. In this paper, we take further advantage of this representation to re-analyse the imposition of the 
simplicity constraints without Immirzi parameter in a GFT context, and thus in the simplicial gravity path 
integral based on BF theory appearing as Feynman amplitude of the theory. We actually solve the above 
two issues, and keep under control the rest of the encoding of simplicial geometry at the quantum level, by 
means of a simple generalization of the GFT formalism itself. This generalization has also the advanta ge o f 
recasting the boundary states of the resulting model explicitly in the form of projected spin networks [2311421] . 
Moreover, we show that the resulting model is unique, in a sense we will specify. The Feynman (spin foam) 
amplitudes turn out to be still a variant of the Barrett-Crane model. We are thus led to re-analyse the 
existing criticisms towards this model. 

While we leave the inclusion of the Immirzi parameter in the formalism to a forthcoming paper [48| . we 
also give a model for the topological sector of Plebanski gravity and compare it to the so-called EPR spin 
foam model. 

We work in a GFT context because it is the most complete setting to study the dynamics of spin networks 
and simplicial structures. But it should be clear that the entire construction could be carried through directly 
at the level of simplicial path integrals (as done in [2l|) or spin foam amplitudes, if one so prefers. 



I. GFT'S IN THE METRIC REPRESENTATION 



In this section, we briefly recall the construction of the metric formulation of GFT's [5|. We start with the 
GFT model for 4d SO (4) BF theory, and then discuss, in the metric framework, the implementation of the 
simplicity constraints Q leading to the (usual versions of the) Barrett-Crane model. 



A. GFT for BF theory 

The GFT for SO(4) BF theory is defined in terms of a field ^1234 := y(<?i, • ■ -54) on four copies of the 
group, satisfying the gauge invariance condition: 

Wi G SO(4), <p{gi,- ■ ■ g A ) = J dhip(hgi,hg2,hg 3 ,hg 4 ) (I) 

The dynamics is governed by the action: 

S = 2 j \-^9i} 4 <Pl234 - gj J [d5i] 1( Vl234 ^4567 ^7389 f962 10 V10 851- 

where [dgi] 4 and [dc/i] 10 denote the Haar measures on four and ten copies of SO(4). The combinatorics of the 
interaction term is that of a 4-simplex (the ten indices i = 1, • • ■ 10 corresponding to its ten triangles), while 
the kinetic terms dictates the gluing rules for 4-simplices along tetrahedra. The quantum theory is defined 
by the perturbative expansion in A of the partition function and, by construction, the Feynman diagrams 
are 2-complexes dual to 4d simplicial complexes: vertices, links and faces (2-cells) of the graph are dual 
to 4-simplices, tetrahedra and triangles of the simplicial complex, respectively. Using harmonic analysis on 
SO (4), the Feynman amplitudes take the form of the Ooguri state sum model [27| . 

The metric formulation [5| of the same model is based on the group Fourier transform of the GFT field. 
The SO(3) group Fourier transform 

dun, 

which maps functions on the group to functions on its Lie 
algebra su(2), straightforwardly extends to functions of (several copies of) SO(4) ~ SU~(2) x SU + (2)/Z 2 ; it 
is invertible on even functions 1 f(g) = /(—<?)■ The group Fourier transform of the field ip is thus the function 



In the sequel we thus assume the further invariance of the Ooguri field under Qi — > —gj in each of the variables. 
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on four copies of the Lie algebra so (4) given by: 

(p( Xl ,..,x A ) = / [dg]V(gi, ..,g 4 )e gi ( Xl ) ■ ■ -e gi (x 4 ), Xi £ so (4) ~ M 6 (2) 



where [dg] 4 is the normalized Haar measure on SO(4) 4 . The plane waves e g : so(4) — > U(l) are the functions 
defined as e g (x) — e lTrx9 , where Tr is trace is the fundamental representation of SO (4). Explicitly, using 
the decomposition of the Lie algebra and group elements x — (x-,x+), g — (g_,g+) into left and right 
components x± 6 su(2), g± £ SU(2), and the parametrization g± = e 9n ' T and x^ —x ■ f of these components 
in terms of the (anti-Hermitian) su(2) generators f— (ti,T2,T3), the plane waves read: 



e g {x) 



where the trace is defined as Trr^ = — 

The space of dual fields inherits by duality a non-trivial (non-commutative) pointwise product from the 
convolution product on the group. It is defined on plane-waves as 

(e g *e g ,)(x):=e gg ,(x), (3) 

extends component-wise to the tensor product of four plane-waves and by linearity to the image of the 
Fourier transform. 

The Lie algebra variables of the dual field (p are interpreted, geometrically, as the bivectors associated to 
a triangle, in each tetrahedron, in the standard discretized version of BF theory. Upon Fourier transform, 
the gauge invariance condition (fTJ) translates into a closure constraint C for the bivectors associated to each 
tetrahedron (field), obtained as the group Fourier transform of the gauge invariance projector V: 



V(p(gi,..,gi) = J dh(p(hgi,..,hgi), -t Ptp = C*tp, (4) 
where C(xi, ■ ■ ■ ,2:4) = £o(Xm=i %n), and Sq is the element x = of the family of functions 8 X defined as: 

S x (y) •= J dhe h i(x)e h (y). (5) 
These functions play the role of Dirac distributions in the non-commutative setting, as 

d 6 y(6 x *f)(y) = f(x) 7 (6) 



where d 6 y is the standard Lebesgue measure on so (4) ~ M 6 . 

The combinatorial structure of the GFT action in the metric representation is the same as in group one, 
while group convolution is replaced by the ^-product. Using the short notation ^1234 := <p(xi, • • • , £4), we 
can write the action as 



S[0] = j [d Xi] ^1234 * <?1234 + A J [d Xi] ^1234 * <?4567 * ^7389 * <?96210 (7) 

where it is understood that *-products relate repeated lower indices as <pi * (pi := ((/>* <p.)(xi), with <p-(x) — 

In this representation, the Feynman amplitudes of the GFT can be computed by sticking together with 
the ^-product the propagator and vertex functions given by: 

„ 4 5 10 

P{x,x') = / [dh] l[6. x M), V(x,x') = / Y[[dh { ] l[(S_ xi * e heh - } )(xf) (8) 
J i=i J £=1 i=i 

The vertex function has twenty Lie algebra variables: two variables xf,xf for each of the ten triangles of 
a 4-simplex, corresponding to the two tetrahedra 1,1' sharing that triangle. The vertex function encodes 
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the identification of the bivectors associated to the same triangle in different tetrahedral frames |5[ , up to 
parallel transport between these frames given by the group elements hihj, 1 ; the sign difference reflects the 
fact, in an oriented 4-simplex, a triangle inherit opposite orientations from the two tetrahedra sharing it. 

The Feynman amplitudes in this representation are simplicial path integrals for BF theory, involving an 
integration over one SO(4) variable hi for each link of the Feynman graph, equivalently of each tetrahedron 
of the dual simplicial complex, which is interpreted as parallel transport between the two 4-simplices sharing 
that tetrahedron; and a single so (4) variable x* for each triangle t of the dual simplicial complex: 

Ag= [Hdhi ntd 6 ^]e iEtlVa:tHt , (9) 

J l e 

dhi is the Haar measure on SO (4) and d e xt is the Lebesgue measure on so (4) ^R 6 . The group element 

df t hi is the ordered product of group elements hi on the links on the boundary of the face (2-cell) 
ft of the graph dual to the triangle t, giving a measure of the discrete curvature associated to it. Such a 
discrete curvature is calculated for a given orientation of ft, and a given 'reference' vertex in its boundary; 
the final amplitude is independent of this choice. The exponential corresponds to the product of plane waves 
Yit Q H t { x t)] its argument Y]. HvXfHt is the unconstraint discrete BF action. 

Importantly, the expression generalizes naturally 5] to the case of open Feynman diagrams (GFT transition 
amplitudes), i.e. simplicial complexes with boundary, to give BF path integrals with fixed B variables on 
the boundary, with appropriate boundary terms (Hoj . 



B. Imposing the simplicity constraints 

This new representation of the Ooguri model for BF theory provides a convenient starting point for 
imposing in a geometrically transparent manner the discrete simplicity constraints that turn BF theory 
into 4d simplicial gravity. The simplicity conditions of Plebanski gravity [28[ (or the gravitational sector 
of Hoist gravity [29(), thus without any Immirzi parameter, ensure that the B-field can be expressed in 
terms of 1-forms fields e/ as B IJ — e IJKL ei A ej. At the discrete level, part of these conditions (in 
their linearized form [22l. l39l are implemented by requiring that, for each tetrahedron in the simplicial 
complex, the four bivectors Bj J associated to its four triangles j = 1, •••4, expressed in the reference 
frame of the tetrahedron, all lie in the hyperplane orthogonal to some normal vector in R 4 . Concretely, 
using the canonical decomposition of bivectors into selfdual and anti-sefdual su(2)-components B ±l Ti with 
2B ±l = \e iik B jk +B 0i , the conditions are imposed by requiring the existence of k £ SU(2)(~ S 3 ) independent 
of j such that: 

Vj, B+ = -kBjk- 1 (10) 

Back to GFT, where the role of the discrete B are played by the field variables Xj in the metric represen- 
tation, the condition (|10[) can be implemented on the GFT field by defining the projector (p \—¥ S k -k@, where 
Sk is the function of four so(4)-variables defined in terms of the non-commutative 6 functions ([5]) as 

4 

S k {x ir -- Xi ) = ,. ; (11) 

3=1 

Such a projector imposes strongly the relations (ITU1) on the field variables, since the non-commutative 5 
functions act as Dirac distribution for the *-product. 
Upon group Fourier transform, one can show that 



Sk * <p = / [dgj] 4 (S k ip)(gjj) e 9l ■ ■ ■ e 94 

^SO(4) 4 



(12) 



G 



where S k projects onto fields on the product of four copies of the homogeneous space S 3 ~ SO(4)/SO(3)fc, 
SO(3)fc being the stabilizer group of k seen as a vector 1 of R ; namely, using the decomposition g = (g~,g + ) 
into self dual and anti-self dual components: 

(S k ip)(g) := / [duj] 4 (p(k 1 u 3 kg-,u j g+). 

The case k — 1 reproduces the standard Barrett-Crane projector [l9|, [2(J. This constraint is then imposed 
on each GFT field, i.e. on each tetrahedron of the simplicial complex, in its associated frame. The task 
of correctly parallel transporting these conditions in the other frame associated to other tetrahedra and 
simplices of the simplicial complex is accomplished by propagator and vertex GFT functions. 

The above amount to, roughly speaking, half of the full simplicity conditions ensuring geometricity of 
the bivectors. The other half, often linked to the secondary constraints arising from the canonical analysis 
[39l [5lj , can be shown to be imposed automatically if the closure condition is imposed on all tetrahedra in 
the simplicial complex (see, e.g. [47]), and if bivectors are correctly parallel transported across frames. This 
condition is imposed on GFT fields already at the level of the action, via the 'closure' projector C given 
by (jU), dual to the gauge invariance projector. Thus, it is imposed on every tetrahedron of the simplicial 
complex. Therefore, we do not impose further conditions on the bivectors in what follows. 

By combining the Barrett- Crane projector S\ with closure, one can build up the field := Si * C * <p, 
defining a geometric tetrahedron. However, one immediately notices that, because the two functions S k and 
C do not ^-commute, the product S k *C does not act as a projector. This fact can be shown to be the 
manifestation of a geometric inconsistency in the way the resulting model imposes the simplicity constraints. 

In fact, given /isSO(4), one has 

(e h * S k )(x) = (S ht>k * e h )(x) (13) 

with h t> k := h + k(h') . This expresses the fact that, after rotation by h, simple bivectors with respect to 
the normal k become simple with respect to the rotated normal h > k. Imposing the simplicity condition 
everywhere with the same normal k simply misses the fact that, under rotation h, normals should be rotated 
as well (which strengthens the relation between the normal variables and the gauge connection) . This means 
that the simplicity constraints are imposed non-covariantly by this procedure. 

But this is exactly this procedure that defines the usual versions of the Barrett-Crane models. In fact, 
combining the interaction term: 

J ^1234 * ^4567 * ^7389 * ^962 10 * ^10 851 (14) 

where \& := Si ★ C ★ 0, with the possible kinetic terms: 

\ J *U34 . \ /(5*£)l234 Or i / (15) 

gives the versions of the Barrett-Crane model derived in [l!|, [2(3] or [2l| respectively. The existence of 
different versions, all differing, in spin foam representations, by edge (tetrahedra) amplitudes only, is then 
understood to be due precisely to the non-covariant imposition of the simplicity constraints (see |50|). 

In this metric representation, the GFT Feynman amplitudes of the constrained theories defined by (|14[) 
and (1151) take the form of simplicial gravity path integrals, with a measure weighted by a constrained BF 



1 Thanks to the isomorphism SO(4) = SU(2) X SU(2)/Z2, each element of SO(4) can be represented as g = (g ,g + ), with 
g ± G SU(2). The map SO (4) -> SU(2), g >->• with kernel SO(3) = SU(2)/Z 2 , realizes the identification of the 

homogeneous space 5 3 ~ SO(4)/SO(3) or set of normal vectors in R 4 , and the group SU(2). Through this identification, 
the rotation by g = (g~ ,g + ) of a normal vector represented by k £ SU(2) is represented by g + k(g~)~ 1 . The stabilizer group 
SO(3)fe of k is thus the set of g = (g~ ,g + ) £ SO (4) such that g~ = k~ 1 g + k, ie of the form g = (k' 1 uk, u), with u £ SU(2). 
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action [5]. A more general version of this calculation is detailed in the next section, so we do not repeat 
it here. The spin foam representation is then obtained by the successive application of the group Fourier 
transform and of the Peter- Weyl decomposition to the same amplitudes, and by the integration of all but 
the representation variables to encode the dynamics of the quantum simplicial geometry. 

One can then analyze the encoding of simplicial geometric conditions in this model directly from the form 
of the GFT action, and test the criticisms mentioned in the introduction. Here we only notice that the 
model couples the bivector variables x across simplices and that the interplay of simplicity constraints and 
parallel transport (following the covariance constraint) leads to further conditions involving both bivector 
and connection variables. It would be tempting to interpret them as a discrete version of the secondary 
constraints of the canonical theory, but this cannot be done in a rigorous way at the present stage. At 
the same time normals to the same tetrahedron, as seen in different 4-simplices, are not correlated, as they 
should on geometric grounds, and this in turn implies a missing geometric condition on connection variables 
hi - the one resulting from the requirement that it correctly transports normal vectors across simplices. 

We postpone a more detailed discussion of the simplicial geometric aspects to the last section. But first, we 
will introduce a modification of the above construction in order to solve the issue of non-covariant imposition 
of the simplicity constraints, and thus the consequent lack of constraints on the discrete connection and of 
correlations among normal variables. This is the goal of the next section. 

II. EXTENDED GFT FORMALISM 

We have seen that some of the problematic issues in the standard formulation of the Barrett-Crane model 
have to do with the way normal vectors to each tetrahedron are treated. A simple generalization of the GFT 
field on which both the Ooguri and Barrett-Crane models are based leads to an easy solution of them. 

A. Tetrahedra normals and gauge covariance 

The needed step is to promote the would-be normal vector to the tetrahedron to a dynamical variable 
k 6 S 3 ~ SU(2), by adding it as an independent argument to the Ooguri field, representing the tetrahedron 1 . 
Thus, in this extended formalism, GFT are defined in terms of fields on so(4) 4 x SU(2): 

<Pk(xir-- > x 4,) ■= <fi(xi,--- ,XA]k) = J [d5]V(5i, • ■ • , 94; k) eg 1 {x 1 )e g2 (x 2 )e g:i {x 3 )e g4 {x4) (16) 

In adopting the new field as our basic dynamical variable in a new formulation of quantum 4d BF theory, 
we still want to impose the local gauge invariance of the simplicial theory, as the invariance under change 
of local frame in each tetrahedron. This time, however, the local rotation should rotate simultaneously the 
bivectors and the normal vectors. In other words, the gauge invariance condition for the Ooguri field is 
now replaced by a gauge covariance of the SO(4) arguments with respect to the normal k, with the overall 
5-argument field being invariant: 

tp(hgi,--- ,hgi;h>k) = (p(gi,--- ,gi;k) Vh £ SO (4) (17) 

where h > k :— /i + /c(/r) _1 denotes the action of the SO(4)-rotation h := (h~,h + ) on the normal vector 2 
k G SU(2) -5 3 . 

To express this invariance in the metric representation, let us introduce the family of functions labelled 
by h € SO(4), given by a product of four SO(4) plane waves e^(x) := e h -(x~)e^(x + ) 3 : 

Ch(xt, ■ ■ -Xi) = e h {xi) ■ ■■e h {x i ) (18) 



1 Obviously, before any geometricity condition is imposed, constraining the bivector variables to be triangle area bivectors and 
relating them to the normal vector, any geometric interpretation is to be understood as purely suggestive. 

2 See footnote on Sec. II Bl 

3 Note that this is the same functions used to impose gauge invariance in the Ooguri BF model. 
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Upon group Fourier transform, the gauge covariance condition (|17[) reads 

(Ch * 0h-ipk){xi, ■ ■ ■ = <?k{x\, ■ ■ ■ V/i G SO (4) (19) 

This condition is implemented by means of the constraint projector: 

C(p k := dhC h *<Ph-i>k (20) 



Two important points must be noted, here. First, the SO(4)-rotation invariance (|17[) of the extended 
field induces an invariance under the stabilizer group SO(3)fc of the normal k, affecting only the first four 
arguments of the GFT field. Upon Peter- Weyl decomposition of the ip^, this means the fields modes are 
labelled by one normal vector fc, four irreducible SO(4) representations (given by pairs of SU(2) spins 
Jj = i = 1 ■ • -4), each of which can be further decomposed into SO(3)fc representations fc»; and a 

four valent SO(3)fc intertwiner contracting these representations. In terms of the dual spin network vertex, 
the labeling thus corresponds to a SO(4) representation in each link, and a normal vector k and a SO(3)fc 
intertwiner at the node. A set of basis functions is thus given by: 



where repeated lower indices are summed over. (g^) are the SU(2) Wigner matrices, (tj) fc< form a basis 
of four- valent SO (3) intertwiners, labelled by an intermediate spin j. The k-dependent coefficients, defined 
in terms of the SO(3) Clebsch-Gordan coefficients C^J^ as: 

C fiJ J k " (k) = £> J * -(k) (22) 

m m pi ' J mm pi mm i y ' v ' 

m 

define a tensor that intertwines the action of SO(3)fc in the representation (g> jf and the action of SO(3) 
in the representation ki. Namely, given = (k^uk, u) € SO(3)&, we have: 

C fi i tk ; (k)D f >_ .{u k )lft. + (u+)=C j l j \ ki (k)D k q i p .(u). (23) 

PT|) corresponds to the vertex structure of the so-called projected spin networks of the covariant approach 
to loop quantum gravity [23| . which thus label the polynomial gauge invariant operators in the extended 
formalism. The convolution of GFT fields to define generic observables, and thus generic boundary spin 
network states at the level of transition amplitudes, will still be defined with respect to the SO(4) connection. 
The need for a generalization of the standard gauge invariance condition on quantum states (Gauss law) in 
the quantization of Plebanski gravity, and in the presence of linearized simplicity constraints, which then 
leads to proj ected spin network states, has been argued from a canonical continuum point of view in [5a ] 
(see also H3 SO, S3). 

The second point is that, despite this generalization, it is easy to check that the following action of the 
invariant (under generalized gauge invariance) field (pk- 

S[Pk] — 7J J ^1234, fc _ 5? y ^ 12 34,fci ^4567,fc 2 ^7389,fc 3 ^962 10,fc 4 VlOSSlfes- (24) 

where we used the shorthand notation y>i234,fc :=<^fc(ffi, ■ • ■ ,54), defines the same amplitudes as the original 
Ooguri model, thus still corresponds at the perturbative level to quantum simplicial BF theory. 



B. Simplicity constraints and geometricity projector 



To obtain a geometric theory in the extended framework, we would like to impose the Plebanski simplicity 
constraints on the model Just like in the previous section, one can implement these constraints by 
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means of the functions Sk labeled by k € SU(2), denned in (fTTj) . where this time the label k is coupled to 
the additional field variable. We thus define a simplicity projector h4 S(p acting on the extended fields 
(pk(xi, ■■■Xi) ■= 0(xi, ■ ■ -xa; k) as: 

(S<p)k := S k *ipk (25) 

Thus the only modification with respect to the GFT formulations of the Barrett-Crane model is that now 
the simplicity projections have to be defined imposing orthogonality of bivectors with respect to the normal 
vector that comes as a new argument of the field. This leads however to a very important difference with 
respect to the previous case. 

In fact, remarkably, this simplicity projector commutes with the action of a rotation, in the sense that, 
given /ieSO(4), 

Ch * (SW^fc = Ch * Sh-ipk * (fh- 1 >k = &k * C/i * tfh-iok — S(Ch * <Ph- 1 >k) 

where ft," 1 > k :— h~k(h + )~ 1 and Ch are the functions of so(4) 4 defined in (ITS)) . The first term correspond to 
the action of a rotation h on the projected field S<p; the last term corresponds to the projection by S of the 
rotated extended field Ch *0h- 1 >k- The simplicity projector acts on the rotated field by imposing simplicity 
with respect to the rotated normal vector. Thus, the equality just states that rotating by h bivectors that 
are simple with respect to a normal k gives bivectors that are simple with respect to the rotated normal 
K 1 >k. 

This commutation property is an important property, because it means that the model constructed in these 
terms will not suffer from the 'lack of covariance' noted in [f| in the various versions of the Barrett-Crane 
model, and that the simplicity constraints are now going to be correctly imposed in the different frames 
across the simplicial complex, obtained as a Feynman diagram of the GFT model in the following section. 

This has a further crucial consequence. Simplicity and gauge covariance are encoded into two projectors 
S and C, defined in (II 1 p and (|20[) . which are now commuting: their product therefore defines a projector 
G = SC = CS, with G 2 — 1. Explicitly, in the metric representation, this projector acts as: 

★ [e h (xi) ■ ■ ■ e h (xi)] * (p(xx, ■ • • x 4 ; /i"fc(ft + ) _1 ) (26) 

Upon group Fourier transform, it acts dually on extended group fields as 

(G^OGfc; k) = f dh ( TT / duA u ,k<i h * */ ,</ * :: ft+fc(/0 _1 ) (27) 

JSO(i) \ jJsU(2) / 

The geometric nature of the tetrahedron corresponding to the GFT field tp is then enforced by the single 
projector G that we can call geometricity projector. It is then imposing this projector on the generalized 
fields that one can construct a GFT (and spin foam) model of 4d quantum gravity based on the (linear) 
Plebanski formulation of classical gravity, free of the mentioned problematic features (from the point of view 
of simplicial geometry) of the usual formulations of the Barrett-Crane model. This is what we do in the 
following section. 

III. A REVISED MODEL FOR 4D QUANTUM GRAVITY WITH NO IMMIRZI PARAMETER 

A. The GFT model and its Feynman amplitudes 

We start from the extended Ooguri model introduced in section [JTJ and use the geometricity projector G 
combining both (extended) closure and simplicity constraints in a single operator. We define a projected 
field \& := G0, representing a geometric (non-commutative) tetrahedron, whose geometry is characterized 
by four area bivectors Xj and one normal vector k. 



(G0)(x 3 ;k) 



dh 



SO(4) 



n 
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The GFT model is defined by the action: 

= \ Jl^Xt} 4 J dfc$^ 34 fc + A y" [d 6 Xi] 10 [dfc a ] 5 *1234,fc„ * *4567,fe 6 * *7389,fc c * *962 10,fc d * $10851, fc e , 

(28) 

where d e Xi is the Lebesgue measure on so (4) ~ R 6 and dfc is the Haar measure on SU(2). It is understood, 
just like in the previous sections, that ★-products relate repeated lower indices as * &i :— * ^f-)(xi), 
with &-(x) = x). 

We can immediately notice four main features of the model so defined: 

• Just as in the Ooguri model and in the standard Barrett-Crane model(s), bivector variables associated 
to the same triangle in different tetrahedra (fields) within the same 4-simplex, as well as across different 
4-simplices, are identified in the frame of the same 4-simplex, to which the above writing refers. When 
expanding the gauge covariance constraint in integral form and writing explicitly the delta functions 
relating the arguments of the fields, one sees that the same bivector variables are now expressed in the 
frames associated to the various tetrahedra, and are related by the discrete connection introduced by 
covariance constraint (see also @, H3, H3, ) . 

• Just as in the extended Ooguri model, normal vectors to tetrahedra are coupled only indirectly in 
each 4-simplex, through their relation with bivector variables imposed by the generalized covariance 
condition and by the simplicity constraints S}~. On the other hand, the normal vectors associated to 
the same tetrahedron in the two different 4-simplices sharing it are strictly identified by the kinetic 
term in the reference frame associated to the tetrahedron (this becomes apparent when making explicit 
the form of the interaction and kinetic functions, and writing the generalized closure constraint in 
integral form). 

• The use of the generalized geometricity projector G imposes all the necessary geometric conditions on 
the bivector variables in each tetrahedron. Note that the closure of the four field variables, hence of the 
bivectors associated to the four triangles of a tetrahedron, is recovered by integration over the normals. 
In fact, such integration makes the insertion of a closure constraint in front of each field ^1234, k in the 
interaction term redundant, since: 



dk5(xi H ^x i )*^(xi,---,x i ;k) = J dk^(xi,- ■ ; x 4 ; k) (29) 

The nature of the operator G as the product of two commuting projectors leads to a unique definition 
of the constrained model. By contrast, the standard Barrett-Crane projector does not commute with 
gauge invariance, giving rise to some ambiguities in the form of edge amplitudes of the spin foam 
model. Here, for example, the use of the interaction term of (|28[) with any of the kinetic terms: 



1234,fc » \ Ji C > £)l234,fc \ J(S k t> <p)wu,k 0r \ J ^1234,fe 



leads to the same amplitudes upon perturbative expansion. At the same time it should be stressed 
that any of the above choices would lead, strictly speaking, to the definition of a different field theory, 
for what concerns the related space of fields, the classical equations of motion, etc. 

The Feynman amplitudes of the model (f2"51) . in the metric representation, can be computed with the 
propagator and vertex functions given by 1 : 

. 4 6 

P(x i ,x' l ;k,k')=6- Xt (x' i )5(k'k- 1 ) V(xix(;k l) k e ,)= / II^II^-*? * e M*-J * S k*')(*f ) (30) 

^ l=\ i=l 



1 This corresponds to the choice: S(ip) = 5/^34^ + 37 / *i234,fc a * *4567,fc,, * $7389, fc c * ^962 lo,k d * *10 851,fe e 
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where Sk(x) := S_ kx - k .i (x + ) and dhg is the Haar measure on SO(4). Note that the vertex does not encode any 
correlation among normal vectors for different tetrahedra. The propagator, instead, presents a fifth strand 
corresponding to the identification of normal vectors in the two 4-simplices sharing the same tetrahedron. 

The calculation of the Feynman amplitudes is similar to the topological case @. In building up the 
diagram, propagator and vertex strands are joined to one another using the ^-product. The amplitudes are 
integrals over the bivector and holonomy variables, and the normals. In terms of the simplicial complex 
dual to the Feynman graph, bivectors x\ are labelled by a couple {triangle t, tetrahedron t} with t C r; 
holonomies h aT are labelled as by a couple {4-simplex ct, tetrahedron t} with rCtr; the normals are labelled 
by tetrahedra as k T . 

Under the integration over the normals and holonomy variables, the amplitude of a closed graph factorizes 
into a product of contributions At[h aT , k T ] for each face of the graph, hence for each triangle t of the dual 
simplicial complex, which take the form a cyclic ^-product: 



A t [K T ,k T ] 



Nt 



3=0 



6^1 



]*^o(^*^*e fcjJ+1 )(^) 



(31) 



The integer j G{0 • • -Nt} labels the tetrahedra sharing the triangle t, or the links of the boundary of the face 
ft dual to t in the Feynman graph. The expression is defined for a given ordering of these tetrahedra, induced 
by a choice of orientation of ft and a reference point in its boundary. hjj + \ := h~Jjh a j + i is the holonomy 
from the tetrahedron j to the tetrahedron j + 1 through the 4-simplex a sharing them. By convention we set 
xn+i '■= xq. The variables x\ G so(4) corresponds to the same bivector associated to the triangle t, expressed 
in the various reference frames of the tetrahedra j. The measure features both the simplicity constraints 
for the bivector at the level of each tetrahedron and the parallel transports of (constrained) bivectors and 
normal vectors across simplices, via the discrete connection. 

The Feynman amplitude can be written as a simplicial path integral, by first integrating, in each face 
contribution At[h aT , k T ], over all bivectors x J t save one Xt '■— xq corresponding to a 'reference' tetrahedron 
frame j — 0; and by a repeat use of the commutation relation eh* Sk — Shi>k * e^, with h t> k = h+kh' 1 , to 
place all the plane waves to the right of the integrand. We obtain the following expression for the graph 
amplitude: 



Ag=j n d/wijd 6 x t /ii dfc -n*^o^>fc,( 



(*0 



* e 



J2 t Trx t H t 



(32) 



where the holonomy h 



Oj 



hoi 



V-ij 



parallel transports the reference tetrahedron to the tetrahedron 



j (with /ioo = 1) &nd H t := ft-oi ' ■ • h,No is the holonomy along the boundary of the face ft of Q dual 
to the triangle t, calculated from the reference tetrahedron frame. The measures dh TlT ,d 6 Xi and dk T arc 
respectively the Haar measure on SO (4), the Lebesgue measure on so (4) ~ K 6 and the Haar measure on 
SU(2). The exponential corresponds to the product of plane waves Y[t eH t( x t)i its argument ^2 t TTx t H t is 
the unconstraint discrete BF action. 

This amplitude takes the form of a (non-commutative) simplicial path integral for a constrained BF theory. 
The measure features the simplicity of the bivector xt of a triangle in the frame of all tetrahedra sharing the 
triangle. Thus, Sh 0j> k,(xt) imposes the simplicity of xt with respect to the hoj >kj, which is the normal to 
the tetrahedron j parallel-transported to the reference frame 0. This is the same as requiring the simplicity 
with respect to kj of the bivector parallel-transported to the frame of j. 

Note that the above amplitude automatically encodes the closure constraint for the bivectors within each 
tetrahedron. This is a direct consequence of (extended) gauge invariance and the integration over the normals. 
Let us nevertheless check it explicitly, by picking a tetrahedron r and considering the variables Xt i: i = 1 • • • 4 
associated to its four boundary triangles, namely, bivectors expressed the reference tetrahedron frame of each 
triangle. The closure constraint says the sum of the four bivectors, once parallel-transported to the frame of 
t, is zero; in the non-commutative formalism, it is imposed using the non-commutative delta- function ([5]). 
Using suitable gauge transformation, we can always suppose that r is the reference tetrahedron for triangles 
U. By appropriately flipping signs Xf t i— > —x^, we can also assume that the orientations of the faces ft t 
chosen to calculate the holonomies induce the same orientation of the common link dual to r. Let us insert 
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the closure constraint <SQ^ xt ; ) := J dg J}j e g (xi) on the left of the simplicity functions in the integrand of 
(|32[) . using the star product. The term depending on x^ in this integrand then becomes: 



i=l 



« 7=0 °'io 



★ e 



d9 II *S> fo) * e*X< Tr ^^ 



where the equality follows from the commutation relation (I13p between plane waves and simplicity functions 
and the definition of star product of plane waves. Now, the dependence upon g in the integrand can be 
removed by the change of variables k T t— > gk T , h aTT i— > h a ^. T g , where a T is the 4-simplex 1 to which both 
t and its successor in the ordered lists j — - ■ ■ N ti of tetrahedra sharing ij, belong. This shows that the 
insertion of the closure constraint is redundant, hence that such constraints are already implemented in the 
simplicial path integral. 

The integrand of (|32[) . as a function of the bivectors, holonomies and normals, is invariant under the 
following gauge transformations, generated by an SO(4) element g T for each tetrahedron: 



h aT ^ h aT g~^, k T H> g+krig'r)' 1 x t >-> g~ r \ t) x t g T ( t ) 



(33) 



where r(t) is the reference tetrahedron for the triangle t and g^ arises from the selfdual/antiselfdual decom- 
position of g T . This corresponds SO (4) rotation of the tetrahedral frames. The integrand is also invariant 
under rotation of the simplex frames acting only on holonomies h ar M- g a h aT (all the other variables being 
expressed in tetrahedral frames). Choosing g T — (<?+, g~ T ) := (1, A;" 1 ) as gauge parameters corresponds to the 
'time gauge', which fixes the normals to the north pole of S 3 . This shows that the integrals over the normals 



k T drop out of the amplitude 
we obtain: 



Making explicit the simplicity functions and denoting hoj := ^(^ch) 



Ag= f JJ d/i CTT []d ( 



X, 



★ e 



(34) 



Note that gauge fixing the normals is only possible for amplitudes of closed graphs; the amplitude of open 
graphs, dual to simplicial complexes with boundary, has an explicit dependence on the normals to the 
boundary tetrahedra. We illustrate this in the case of the simplex boundary state in Sec. IIII CI 



To summarize, we have written the GFT Feynman amplitudes as simplicial path integrals with the explicit 
form of a constrained BF theory, with a clear form of (linear) simplicity constraints imposed in each tetra- 
hedral frame, with respect to a given normal variable and correctly parallel transported across frames by 
means of the gauge connection. The integration over the same connection imposes the closure constraint in 
every tetrahedron; this in turn implies the remaining (volume) simplicity constraints [22, 47]. The amplitude 
also encodes constraints on the connection {h aT }, following from simplicity, giving a specific measure on the 
space of connections. 



B. Spin foam representation of the amplitudes 

The same amplitudes can be written in spin foam representation, exploiting the exact duality between 
the metric and spin representation of the GFT [5]. This re- writing is obtained by inverse group Fourier 
transform - giving the pure lattice gauge theory formulation of the model involving only group elements [J] ) 
- and the Peter- Weyl decomposition of the amplitudes, leaving only a sum over representation of pure spin 
foam amplitudes. 

For the present model the calculation goes as follows. We first notice that, under the integration of the 
group elements h aT in (|34p . each integral over Xt takes the form of an inverse SO (4) Fourier transform 



1 cr T is also the 4-simplex the dual to the 'target' vertex of the oriented link dual to r. 
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formula. Indeed, for a given set of holonomy variables {/i CTT }, the functions "A^Iq ^-h 0jX ih^X x t) m side the 
brackets, labelled by t, define functions of su(2), image by non-commutative Fourier transform of some 
functions O t on SO(4): 

*f=o Z-H 0jX - t h- ) to) = / dgO t (g)e g (x) 



SO(4) 



The inverse Fourier transform formula (see for e.g @) gives: 



(35) 



A quick comparison between this formula and the integrand of (|34[) leads to the following simple expression 
of the graph amplitude in terms of the group functions Ot ■ 



Ag = J Y[dh Ta Y[O t (H^) 



(36) 



where we recall that H t :— hoi ■ ■ ■ h^o is the holonomy around the triangle t. 

Let us determine the functions Ot- The plane- wave expansion of the non-commutative delta- functions 
gives: 

*U 5 r hojX -- K) {xt) = [dgj n d ^ e 9 nf =1 ^fto/^-) e 9 nf =1 < 93 ^; l(x+) 

3=1 

= J n d -9±/ II d * s(9+9?>g. nf =a h- 0j9j h- f e 9+ nf =1 < J -*-^; 1 (*+) ( 37 ) 



where the integrations are over SU(2). We omitted the t of N t to simplify the notations. After the following 
change of variables: 



JV 



9- >-> 9- II h' 0j gjh- -} , 

3=1 



AT 

3+ ^ 5+ IJ KjgjKj 1 

3=1 



(38) 



the expression (|37|) takes the explicit form of a SO(4) Fourier transform Ot(x) — J dgOt(g)e g (x), where Ot(g) 
is given by: 



N 

Ot(g)= I Y[d 9j S( 

3=1 



3 = 1 



AT 



3=1 



(39) 



The evaluation of this function for g — H t 1 gives: 



f N i 

J 3=1 



(40) 



Next, we use the Plancherel decomposition of the SU(2) delta function in terms of the characters x J (ff) = 
TrJ J (q) in the SU(2)-representations labelled by J G |N, which also labels the simple SO (4) representation 
(J, J) Ullil, to write: 



OtiHt 1 ) = d J / II d #J' X J ([Ki9i h i2 ■ ■ ■ gNh^}' 1 h~ 0igi h~ 12 ■ ■ ■ g N h~ N0 ) 
j j=i 



(41) 
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where dj — 2J+1. The integration over each group variables gj in (|40p. which appears exactly two times in the 
argument of each character x J \ is performed using N times the orthogonality relation of the representation 
matrices 



In fact, if we denote Xo '■= X' J (^01(^01) 1 an dj for if = 1 • • • iV: 

Xk '■= X J ([Ki9iht 2 ■ ■■9Kh^ K+1 ] 1 h' 01 gih' 12 ■ ■ ■ g K h~ KK+1 ), 
(where by convention N + 1 =0), one can show, using orthogonality, the following recursion relations: 

/ A9kxJr = jjX^WC^JC+i)" 1 )^-! ( 42 ) 

We conclude by iteration: 

OtiHf 1 ) = E "wvTT II x'CWW 1 ) (43) 

J a J 3=0 

where iV + 1 is the number of tetrahedra sharing the triangle t. 

Using the form (|36| of the Feynman amplitude (|34|) . and the above expression for the functions Ot, we 
finally obtain the spin foam amplitude: 

^ = e n d i n tt^7 n^ 10 j ^ ^ 

{Jt} t T llte9r a J t £7 

The sum is over the SU(2) spins Jt labeling the triangles of the simplicial complex dual to Q and the 
products are over the triangles t, the tetrahedra r and the 4-simplices a. The 4-simplex weight { 10 J} a is 
the Barrett-Crane lOj-symbol [44j |: 

{ioj} CT =/n d ^ n x^HKiK')- 1 ) (45) 

v l = (vv') 

In this expression, v label the five vertices and I the ten links of the graph dual to the boundary of the 
4-simplex: a vertex is dual to a boundary tetrahedron, a link I — (vv') is dual to a triangle t(l) sharing two 
tetrahedra. The measure is over SO(4) elements h v := (/i+, h~ v ). 

The model is thus a variant of the Barrett-Crane model with a specific edge amplitude, different from that 
of other versions of the Barrett-Crane model [5(| (notice in particular the absence of the norm of the BC 
intertwiner from the amplitudes), although identical to that obtained in [2lj |. whose construction is indeed 
close to ours regarding both its motivations and techniques employed. 

From both the construction and the detail of the calculation, it should be clear that, while the appearance 
of the lOj-symbol is the direct result of the form of the simplicity constraints chosen (corresponding to the 
pure gravity sector of the Hoist action, or to the Plebanski formulation of gravity, with linear constraints), 
most of the information concerning the details of the simplicial geometry behind the model, in particular the 
extended closure condition and associated Lorentz covariance of quantum states, as well as the correlations 
of normal vectors across simplices, is now hidden in the edge amplitudes. The reason for this non-transparent 
form is obviously the fact that all the geometric variables in which this information is manifest, i.e. bivectors, 
gauge connection and normal vectors, have been integrated out. 



C. Simplex boundary state 



Taking the GFT Feynman amplitudes as the definition of a lattice model, we would like to discuss here 
the simplex boundary state. The vertex function gives the amplitude for a single 4-simplex (represented 
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by the vertex graph) with fixed boundary data {x^,xl } and {k T }: the pair xj,xj both corresponds to 
bivectors associated the triangle t, expressed in the frame of the two tetrahedra sharing it. Given ten of 
such pairs, the vertex function encodes the relations between the bivectors of the same triangle expressed in 
different tetrahedral frames: 

fU dh r II {^*SK*S KT , >K ,*e Kr ,){xt), (46) 

r *=<tt'> 

The product is over the triangles of the 4-simplex, each of which shared by an (ordered) pair of tetrahedra. 
h TT > = h~^h T i is the parallel transport between r and r' and Sk[x) '■= S_ kx - k -i (x + ) are the simplicity functions 
written in terms of the non-commutative delta functions ([3]). 

An important remark is in order, here. Without the simplicity functions, the bivectors x\, x\ would be 
identified up to (a sign and) parallel transport h TT / from rtor', as they should be classically. However, due to 
the specific properties of the star product, the presence of the simplicity functions relaxes this identification. 
Using the plane waves expansion of S k and the definition of the star product on plane waves, one can indeed 
show that: 



(5. xI * S K * Sh TT ,>k T , *e KT ,)(xJ ) = j &u T t &u T t (5- xl *S kT *S KtI>k , * e a[ ^ u[ 0(x[ ) (47) 

where the integration of the right hand side is over the product SO(3) T x SO(3) T < of stabilizer groups of the 
normals k T and k T >. The fact that the amplitude does not force the complete identification up to parallel 
transport of the bivectors accross neighboring tetrahedra, also noticed in [21] and argued to be the main 
shortcoming of the Barrett-Crane model, appears, in our formalism, to be due to the star product structure, 
and thus the quantization map chosen for the constraints. As we argue also in the discussion of Sec. |V1 
such relaxation of parallel transport could be purely quantum effects, devoted to disappear in a suitable 
semi-classical limit. In fact such a requirement serves as a guide for the very definition of such a limit, which 
should reach a regime in which the star product is approximated by a commutative product. A scale A > 
can be introduced in the amplitude via a modification of the group Fourier transform, where SU(2) group 

elements are parametrized as g = e x ^' T withpgR 3 and \p\ < 1/A, plane waves are replaced by e^(a;) :=e^ jTrx9 , 
and the star product preserves the parameter (l4| : 



In the limit A — > 0, corresponding to a regime of low curvature, the star product can be approximated by 
the standard pointwise product of functions on so(4)~R 6 . In this parametrized framework, the parameter 
A appears not only in front of the action in the path integral, making it suitable to define the semi-classical 
regime, but also as a deformation parameter of the product used for the definition and the imposition of 
the constraints; this contributes to further quantum corrections, of which, the above relaxation of parallel 
transport may be a manifestation. This possibility will be studied in more detail elsewhere. 

Let us also point out that, if one defines, in the parametrized framework, the su(2) non-commutative delta 
function as 6 — jg f dge^, then they formally become true delta function on su(2) ~R 3 in the limit A — > 0. 
With simplicity functions defined as before Sk(x) = 8_ kx - k -i {x + ), this means that, in this limit, the constraints 
in (|46[) impose, for each triangle t, the following condition on the bivector xt ■= x T t and the connection: 

x^ — k-j-Xj.k"r f — ^ k-p' J x ^ [^tt' t> j 

In words, the normals of the tetrahedra sharing t, expressed in the frame of r, differ from an element in the 
stabilizer group of the x~ t , namely [h TT rt>k T i]k^} 6 U(l) a ,+ This expresses the geometrical fact that the two 
normals belong the plane co-orthogonal to bivector x t (namely orthogonal to his Hodge dual, which defines 
the plane spanned by the triangle). In time gauge k = 1, this becomes a condition on the boost part of the 
parallel transports h TT i between adjacent tetrahedra. 



The integration of the amplitude (|46|) over half of the bivectors variables, namely, for each triangle t = 
(t, t'), the variable x\ gives a function over the normals and the ten remaining variables Xt := xj* , which 
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are bivectors of the triangle t in its reference tetrahedron r t . The resulting amplitude takes the form of a 
discretized (and non-commutative) path integral for constrained BF theory, for a single 4-simplex, with fixed 
boundary bivectors: 

A(x t ; kr)= f J] [dg t e^ Xt9t * S krt (x t )] * / J] dh T [] S Kt , >k , (x t ) S(h T gtK)) (48) 

J t r t 

In this expression, we made explicit the part of the constraint implementing the linear simplicity of the bivec- 
tors in their reference tetrahedron, and the part (also dependent on the bivectors) entering the definition of 
the measure over the connections. Note that the closure constraint for the bivectors is only implemented after 
integration over the normals. That a proper spin foam model should have this form, with the corresponding 
constraints modifying the measure over the connections, has been advocated in [39| . 

IV. A MODEL FOR THE TOPOLOGICAL SECTOR OF PLEBANSKI GRAVITY 

In this section, we sketch the construction of a GFT model for the so-called topological sector of Plebanski 
gravity as a constrained BF theory, where the constraints on the B-field impose that it can be expressed in 
terms of 1-forms B IJ = e 1 A e J rather than e IJKL ei A e,j. At the discrete level, their are implemented by 
requiring the same condition as in (|TU)) . but for the Hodge dual to B rather than B itself. In terms of B, the 
condition differs from ([TU)) by a sign: there exists a k € SU(2), such that Vj, 2?+ = kB^k' 1 . 

This condition is obtained from the BF GFT in its extended formulation, by using an operator, acting on 
the extended field (in metric variables) as (pk >— > (S top ip)k := S k ° p * tp/., where S k ° p are the functions of four 
so (4) variables, labelled by a SU(2) element fc, defined in terms of the non-commutative <5 functions ([5]) as 

S t k op (x 1 ,---x i ) = f[S kx - k - 1 (x+) (49) 

2=1 

Upon group Fourier transform, 



Sl op *0 k = [dgj] 4 {Sk P ipk){9j ) e Sl • • • e 94 

JSO(4) 4 

this simplicity operator acts as: 

0Sj>*)GO = II / t d ^] 4 <Pk(k- l ufkg-,u jg +). (50) 

where we used the decomposition g = (g~,g + ) into self dual and anti-self dual components. 

Note that, unlike its analogue in the gravitational sector, the simplicity operator does not define a projec- 
tor 1 . One can however check that it commutes with extended gauge invariance projector C in (TTT)) . This is 
because, just like their analogues Sk, the simplicity functions S l k ° p satisfy the property that, given h G SO (4), 



*^° P = ^ fc *e, (51) 



This property allows to impose the simplicity constraint covariantly. 

A model can be defined by the action having the same form as ([2"5]). but now defined in terms of fields 
# = G top (p, with G top = S top C = CS top . Just as in the gravitational sector, the Feynman amplitudes take 



1 There are ways to impose such constraints via a projector, for e.g by acting on the field with S k ° p by right star multiplication 
on the components x~ , left on the components x + \ we will not investigate them in this paper. 
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the form of simplicial path integrals for a constrained BF theory: 



Ag = I Yl dh ar ]Jd ( 



Xi 



n 



"j=0 u h aj x t h 



* e 



J2 t Trx t H t 



(52) 



where the notations are the same as in the previous section: for each triangle t, the integers j = • • • Nt 
labels the tetrahedra sharing t, hoj is the holonomy between the tetrahedra and j and h := h + (h')' 1 . Ag 
is the generic amplitude of closed graph expressed in time gauge, as (extended) gauge invariance allows to 
remove the dependency upon the normals. 

The spin foam model can be calculated directly from its simplicial integral form, following the same route 
as in Sec. IIII B[ We will not describe in detail here the details of the model. Let us however have a look 
at the action of the simplicity operator on the Peter- Weyl components on the fields. Gauge invariants fields 
ifk expand into projected spinnetwork vertex functions (|21[) : for k = 1, 
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(53) 



We deduce from the orthogonality relations of the Wigner matrices 

duD^u^D^u) = -j+6 j+t:r 6 m , q 6 ni p 

that simplicity (|50p projects onto simple representations j£ = , just like in the gravitational sector, and 

exchanges the indices m and v! i in C^Jj^ in the formula (|53|) . Given that the Clebsh Gordan coefficients 
satisfy: 

njjM — (_-\ \2o+kfij,j,k 

we conclude that S top projects onto simple SO(4) representation J = (j, j), does not impose any restriction on 
the expansion of (j, j) into SU(2) irreducible k — 0, • • • 2j, and acts on each component ( J, k) by multiplication 
by the phase (— l) 2j+fe . The expansion into the SU(2) representations k is what defines the fusion coefficients 
charaterizing the new models [22|], and in particular the EPR model. In the EPR model, the expansion is 
restricted to spin 2k~j. The model obtained here is manifestly different; the action of the simplicity function 
rather suggests a form analogue to the Ooguri model, as a sum over SO (4) representations of products of 
(15j) symbols, where the representations are restricted to be simple. 



Thus, by following the same strategy as in the previous sections, a model for the topological sector of 
Plebanski gravity can be proposed in the non-commutative metric formulation of GFT, where the classical 
constraints are imposed via non-commutatice delta functions. Just like in the gravity case, the amplitudes 
are simplicial path integral with a clear geometric content. The resulting spin foam model is manifestly 
distinct, at least in the full quantum regime, from the one we would obtain by imposing the constraints on 
coherent states parameters |22j |. 



V. DISCUSSION 



In this section we summarize our results and re-examine in light of them the arguments against the 
Barrett-Crane model, by which, from now on, we mean the amplitudes associated to a given triangulation 
as obtained within the generalized GFT formalism presented above, whether written in the simplicial path 
integral or in the spin foam representation. 

First of all, let us recapitulate the procedure we adopted, and the ingredients we put in. 
We started from a GFT formulation of quantum BF theory in 4 dimensions. The GFT field represents the 
quantization of the data associated to a single tetrahedron in a simplicial decomposition of spacetime (to 
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be generated by the GFT perturbative expansion) in BF theory. Its arguments represent the corresponding 
simplicial BF variables. We extended the standard formalism to one including, for each tetrahedron, an 
additional variable valued in S 3 , in preparation for the imposition of the geometric constraints that should 
turn this model into one for 4d gravity. We then imposed onto such extended GFT field an extended 
covariance condition affecting both the usual field arguments and the additional S 3 variable. This condition 
represents the gauge invariance constraint of (simplicial) BF theory, and introduces, when written in the 
form of an integral operator acting on a generic GFT field, the bulk gauge connection of the theory, encoding 
parallel transport among tetrahedral and simplex frames. In a geometric theory, this should also encode the 
closure constraint on the faces of the same tetrahedron. 

After appropriate constraints, in fact, the GFT field should turn into a quantum description of a (quantum) 
geometric tetrahedron. The needed (simplicity) constraints can be split (both in the discrete and in the 
continuum setting) into two sets: diagonal and cross-diagonal simplicity, acting at the level of individual 
tetrahedra, and here imposed in their combined linear form [221.1471] . and 'volume'simplicity constraints. The 
variables on which these constraints have to be imposed are exactly the Lie algebra elements appearing as 
arguments of the GFT fields in the metric representation, as confirmed by the GFT amplitudes written in 
the simplicial path integral representation. We then impose the first set of simplicity constraints, as well 
as the extended covariance (thus closure) as restrictions on the GFT field. Such restrictions then show up 
in the form of (non-commutative) delta function insertions in the simplicial BF path integral, at the level 
of GFT Feynman amplitudes, for all tetrahedra in the associated triangulation, and in the frames of the 
tetrahedra, thus result in a covariant imposition of the constraints. Moreover, we know [22| . |47| that the 
volume simplicity constraints follow automatically when parallel transport conditions among tetrahedral 
frames, simplicity and closure conditions hold. In turn, when all of these conditions hold, an assuming 
non-degeneracy, the set of bivcctors define a unique geometry for the simplicial complex, and thus can be 
inverted for the set of edge lengths. No classical geometric condition is thus missing in the construction. 

The next step, however, is to quantize these conditions. As mentioned, this process is necessarily ambiguous 
and a choice of quantization map must be made. And as mentioned, our procedure relies on the non- 
commutative Fourier transform and on the resulting star product used to compose the building blocks of the 
simplicial path integral, including the simplicity constraints and the parallel transport conditions. Thus we 
have implicitly used, in the path integral setting, the Duflo quantization map. Once more, support for this 
choice, beside the more abstractly mathematical one, comes from the results of its application to pure BF 
theory |5|. [3^I as well as to simpler systems (33j (and, more recently, to canonical 3d gravity with cosmological 
constant |45j). 

In this quantization, the closure condition of bivectors associated to the faces of each given tetrahedron 
is not manifest, in contrast to the standard formulation of BF theory, due to the extended form of gauge 
covariance. One may then want to include it explicitly, inserting a non-commutative delta functions at the 
level of the Feynman amplitudes of the model. However, as noticed also in 30], and remarked above, it is 
immediate to see that such insertion would be redundant and it could be re-absorbed in the integration over 
the discrete connection. The closure constraint is indeed already implemented by this integration. 

As we have shown in the previous sections, this quantization procedure also implied a characteristic 
interrelation between simplicity constraints and parallel transport, which leads to the covariant imposition of 
the first, resulting restrictions on the connection, and a consequent relaxation of parallel transport conditions. 
This fact is not affected by the presence or value of the Immirzi parameter, but once more is a necessary 
consequence of the quantization chosen. Actually, it seems to go beyond it and be a feature of the other 
spin foam models as well. The result is that the configurations summed over in the simplicial path integral, 
and thus in the spin foam model, are not strictu sensu classical geometries, due to unavoidable quantum 
corrections (it can also be understood to be a consequence of the specific path integral measure that our 
GFT procedure defines). 

The end result, in this case, is a simplicial path integral with non-commutative variables, and in spin 
variables, a specific version of the Barrett-Crane spin foam amplitudes. 

Let us then reconsider the issues concerning the quantum simplicial geometry behind the BC model, includ- 
ing the one mentioned above, as made manifest in the metric/simplicial gravity path integral representation. 
These amounted to: 

a) 4-simplices speak only through face representations, i.e. triangle areas [25j : 
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b) bivectors associated to same triangle in different simplices are not identified (after parallel transport) 

c) normal vectors to the same tetahedron seen in different 4-simplices are uncorrelated [5|, |43[ ; 

d) the simplicity constraints are imposed in a non-covariant fashion [f|; 

e) because of this non-covariance, there are missing constraints over the connection variables Q . 

The extended GFT formalism (with linear constraints and no Immirzi parameter) that we presented were 
meant exactly to solve the last three problems above, and they do. 

The other points refer to the coupling of simplices and to the coupling of tetrahedra within simplices. The 
non-commutative metric formulation shows that the bivector variables associated to the ten triangles in each 
4-simplices are correctly identified across different 4-simplices sharing the same tetrahedron (and thus the 
same triangles), thus addressing the first issue. The additional coupling of normal vectors, in turn covariantly 
related to the bivectors by the simplicity constraints and by the generalised closure relation, ensures that the 
correlations among simplices resulting from bivector identifications are not undone by the lack of correlations 
of normals or by missing conditions on the discrete connection. Obviously, if one expands the transition 
amplitudes by Peter- Weyl into (simple) SO (4) representations and then integrates out all variables except the 
representations themselves (interpreted as quantum numbers for the areas of the triangles), the amplitudes 
associated to individual 4-simplices will share only such remaining variables with one another. This simply 
means that the whole simplicial geometry behind the model, and the exact correlations among geometric 
variables (bivectors and normals) in the same tetrahedron as seen by different (neighboring) 4-simplices have 
been encoded in a not-at-all-transparent way into the lower dimensional amplitudes. 

As we illustrated in Sec. IIII C[ point b), argued in [2lJ, to be a serious shortcoming of the Barrett-Crane 
model, has to do exactly with the relaxation of parallel transport resulting from the interplay between 
simplicity constraints and gauge covariance, due to the properties of the non-commutative star product. 
In fact, using the decomposition of the simplicity constraints into plane waves, the star product allows to 
recombine them with the plane wave corresponding to the unconstraint BF action, into a single effective 
action for the model. This is the action used also in plj . Now, as observed in this work, the equations of 
motion for action do not force the complete identification up to parallel transport of the bivectors across 
neighboring tetrahedra. Is this a problem or a feature of the model, i.e of the quantum amplitudes? We argue 
that it is a necessary feature - which will be present also in any extension of the same construction involving 
the Immirzi parameter-, an indication of what quantum geometry may be, and not necessarily the sign of a 
problem with it. Our view is that the effective action obtained using the star product should be interpreted 
as an effective quantum-corrected action indicating properties of the quantum configurations summed over 
in the path integral, which indeed should not be expected to be classical simplicial geometries, in general. 
As a consequence the variations of the same action should not necessarily be interpreted as encoding the 
classical geometry behind the model, and their failure to reproduce classical simplicial geometry should not 
necessarily be interpreted as indicating that geometric conditions arc missing among the ingredients leading 
to it. It just indicates how these necessary ingredients are affected by the quantization procedure. While we 
leave a detailed analysis of the semi-classical expansion of the amplitudes for future work, we can already 
argue that such expansion will involve approximating the star product with a commutative product, and thus 
a commutative limit of the equations of motion resulting from the effective quantum action appearing in the 
path integral. This is exactly what happens in the simpler case considered in [33j . In this approximation, 
the issue with parallel transport of bivectors disappears and classical simplicial geometry would then be 
recovered. 

In light of the above, the criticisms of (the first reference of) [22|], saying that bivectors are not identified 
across simplices and that they are in the new models, therefore, can also be reconsidered. At a closer 
examination one notices that what is not identified across different simplices in the BC amplitudes are the 
coherent state parameters associated to each triangle in the different 4-simplices, in quantum BF theory, 
on which a quantum version of the simplicity constraints is then imposed strongly (i.e. by means of delta 
functions) in the FK model (and in its extension to finite 7). The point is, however, that the coherent 
state parameters can be identified with the continuous bivector variables of BF theory, only in a semi- 
classical sense i.e in the sense of mean values 22] or in the asymptotic regime of the spin labels 
[22I |53| . This is confirmed by explicit calculations of the group Fourier transform of Wigner representation 
matrices [H, HU, HU , in the coherent state basis, relating the coherent state parameters with the classical 
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Lie algebra variables (for each SU(2) component). This gives, for large spins j: -Di-(x) ~ S(x — jn). It is 
directly the continuous bivector (Lie algebra) variables, on the other hand, that appear in our simplicial path 
integral representation of the amplitudes, and on which the simplicity constraints are imposed, by means 
of (non-commutative) delta functions, as appropriate in a path integral quantization [B|, [24], [3(|. As for 
the identification of bivectors across tetrahedra within a given simplex, the same remark applies to the FK 
model (the other available one in absence of the Immirzi parameter): no strict identification up to parallel 
transport is imposed on the Lie algebra B variables, and, as shown in (2l| . even the identification of coherent 
state parameters is imposed only in a semiclassical limit. 



We now move on to discuss the other arguments that have been put forward against the Barrett-Crane 
model. We do so for the sake of completeness, only. We have not dealt directly with them in this paper, nor 
worked in the context in which these concerns were raised. We limit ourselves, therefore, to reconsider the 
evidence suggesting these concerns, in light of what we learned from our results. 



As a preliminary step, we note a very general criticism raised |24j,[30J against the BC model, which concerns 
other spin foam models as well. While the aim is to obtain a quantization of a classically constrained theory, 
i.e. the Plebanski formulation of gravity as a constrained BF theory, it results from a procedure that 
amounts to first quantizing the classical configurations of BF theory to get its quantum states (SO (4) spin 
networks) and then constraining such quantum states (usually in representation space). However, our metric 
representation of GFTs and spin foam models, providing their dual simplicial path integral representation, 
makes possible to follow a more standard path integral procedure and avoid some ambiguities, and to keep 
under control which of the constraints are imposed and how, even though it still (necessarily) depends on a 
choice of quantization map, as we discussed 1 . 

A popular criticism of the BC model is that it results from imposing the simplicity constraints 'too strongly' 
at the quantum level (22|. Indeed, when the constraints to be imposed are second class, as it is the case in the 
presence of Immirzi parameter or for quadratic simplicity constraints in general, they have to be imposed 
weakly in the quantum theory. On this, we notice two points. First of all, the model presented above 
imposes the simplicity constraints at the level of simplicial path integrals, by means of delta functions, as it 
is appropriate to do both for first and second class constraints. The nature of the constraints is manifested 
in their relation with other constraints and in the resulting further modification of the quantum measure. 
Second, while quadratic simplicity constraints are indeed second class, in the pure gravitational sector of 



Hoist gravity or in the Plebanski theory with linear constraints, the constraints are actually first class 41], so 
that they have to be imposed strongly. The only kind of 'weakening'of constraints that is consistent with the 
classical theory in our case is the one following from the non-commutative nature of the bivector variables 
themselves (whence the dependence on a quantization map 2 ), which is taken into account by the ^-product 
and associated group Fourier transform. 

A related criticism [39|, H(| is that the BC model does not seem impose the secondary second class con- 
straints that arise in the canonical analysis of the continuum theory |55l . l5rjj (again when using quadratic 
simplicity constraints). This criticism applies as well as to all the new spin foam models [22]- It is difficult 
to test any spin foam model in this respect, because none of them is derived following a canonical quan- 
tization procedure. Moreover, being based on a simplicial discretization, the standard canonical reasoning 
is difficult to apply and discrete counterparts of continuum secondary constraints are very difficult to de- 
fine. A variant of this criticism adapted to the simplicial setting comes from the beautiful analysis of [5lj . 
and from the standard discretization of the simplicity constraints of the classical theory, whether linear or 
quadratic [13, [H, E3] ■ These analyses show that in the canonical setting it is necessary to impose all sim- 
plicity constraints, including so-called volume or edge simplicity constraints, on top of the ones we discussed 



1 The same motivation is shared by the work of 13011 . where however the non-commutativity of bivector variables, following 
from their conjugate nature to the non-abelian discrete connection and consistent with the canonical phase space of both BF 
theory and loop quantum gravity, is neglected. 

2 One may also wonder why this quantization prescription is already present in the classical GFT action; this however is simply 
a confirmation of its interpretation as a "third quantization" formalism 15131 . 
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explicitly, and also so-called 'gluing constraints'which are analogue of secondary continuum constraints, in 
order to ensure geometricity of the classical configurations appearing in the model. This is indeed a reason 
for concern, and a detailed analysis is certainly needed in the context of the model we presented, and of 
the other available spin foam models. However, we emphasize that the non-commutative path integral for- 
mulation of spin foam models, to be developed also for the newer models, allows to study in a very direct 
way which constraints one is imposing in the definition of the spin foam model. Second, it is known [22T 
that, when imposed at the covariant classical level in all tetrahedra of the triangulation, the combination of 
diagonal, off-diagonal and closure constraints is enough to ensure full geometricity of the degrees of freedom 
included in the model, modulo degenerate configurations, if gauge covariance (parallel transport) is also 
imposed. Third, we have noticed how our covariant imposition of simplicity constraints involving the extra 
data encoded in the normals and the generalized form of the closure conditions imply further conditions on 
the normal vectors and on the discrete connection variables. The extra conditions affect the propagation of 
simplicity constraints from frame to frame across the triangulation, and thus their correct imposition in all 
tetrahedra, including those lying on different hypersurfaces. Thus it could be that such extra conditions are 
equivalent to the canonical secondary constraints. However, further work is needed to confirm or refute this 
possibility. 

Several criticisms concerned the boundary states of the BC model. Once more, more work is needed to 
define properly the boundary Hilbcrt space of the theory we presented. Still, our results offer some insight. 
It has been argued [22|, 5(| that states BC model do not have enough parameters to correspond to a 3- 
geometry. In other words, there would be a too-restricted set of commuting observables and in particular no 
label for the 3-volume (no intertwiner degree of freedom). In this respect we can notice that the boundary 
states of the BC model, as we have defined it, are given by graphs labelled by SO (4) representations (or 
Lie algebra elements (fluxes) [3l| or group elements) on their links, subject to simplicity constraints, and 
one normal vector k for each node, subject to our generalized gauge invariance condition. They are indeed 
projected spin networks as we have noticed [23l. I39l.l40l.l42j. A heuristic interpretation of such state is as the 
discrete analogue of canonical wave functions on the extended gravity phase space, including lapse and shift 
vector: ty(h,N,N l ). The standard states of the BC model, labelled only by simple representations and by 
the BC intertwiner, arise after gauge invariance is imposed and normals are integrated out. However, this 
integration over normals cannot be done freely (i.e. on each vertex separately, in the boundary spin network 
graph) because of the requirement of covariance, and because normals are correlated across tetrahedra by 
means of the discrete connection (also subject to constraints, as we have seen). It would rather seem that 
the states of the model, now including the normals, possess all the geometric information to be expected, 
while it can of course still be true that particular linear combinations of states (as resulting for example from 
arbitrary integrations of normal degrees of freedom) do not correspond to well-defined boundary data. 

Clearly, this state space does not match the one of loop quantum gravity. This is based on SU(2) spin 
networks, that are labelled by an arbitrary SU(2) intertwiner, and that, while embeddable in a covariant 
way in SO (4) using indeed the projected spin network formalism [39|, H2, 5f|, are based on the Ashtekar 
connection that requires for its definition the introduction of the Immirzi parameter [l|, [22[ ■ Indeed, the 
stated failure of the BC model to reproduce the tensorial structure of the (lattice) graviton propagator [46[, 
at a closer examination, could be interpreted simply as a mismatch between the data encoded in LQG spin 
networks and used to define both the background boundary states and the same graviton observables, and 
the boundary states (and consequent dynamical data) of the BC model, that, as we have shown, should be 
taken to be the above projected spin networks. However, no such matching can be expected for a theory 
that is independent of the Immirzi parameter from the start. 

Finally, the inclusion of the Immirzi parameter is the most interesting feature of the new spin foam models 
[22| , A generalization of the non-commutative metric formalism presented here to include the Immirzi 
parameter is possible and will be presented elsewhere [48j . together with the resulting new model. Still, two 
general points must be noted. First, one could expect any new model constructed in this way to reproduce 
(one of) the new spin foam models only in a semi-classical approximation, having in mind both our results 
on the topological sector of the theory and the construction of the new models based on group coherent 
states. Second, at present our main reason to want the Immirzi parameter in a quantum gravity theory is 
its necessity in the LQG context. It would be good to have further confirmations, either at the classical or 
quantum level, of its importance to capture the correct physics of a quantum spacetime. 
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Last, we note that an important open issue is the role that degenerate geometries play in all spin foam 
models, including the one we have presented. Not only they affect drastically the imposition of simplicity 
constraints, and thus the geometric nature of the configurations included in them [l2l.[5l|. but they may end 
up dominating the quantum dynamics on entropic grounds [26j]. It is also crucial to stress, however, as done 
in [46( | that this issue has to be studied in the context of the computation of physical quantities and that 
the insertion of appropriate observables or boundary states can change dramatically the respective weight 
of proper and degenerate quantum geometries. 

Conclusions 

We have used the non-commutative metric formulation of group field theories and spin foam models to 
define a model of 4-dimensional quantum gravity as a constrained BF theory, without Immirzi parameter. 
This involved a generalization of the usual GFT framework to include both Lie algebra or group elements, 
associated to triangles in the triangulation, and normal vectors associated to tetrahedra of the same. The 
generalization led naturally to projected spin network states and the associated covariance under Lorentz 
transformations. The model is uniquely defined thanks to the projector nature of the generalized geometricity 
operator we introduced at the GFT level, encoding both simplicity constraints and (generalized) gauge 
covariance. The resulting model, for which we exhibited both a complete simplicial path integral expression 
and a spin foam representation, turns out to be a variant of the Barrett-Crane model, characterized by 
specific lower-dimensional amplitudes. We also presented a similar construction for the topological sector of 
the Hoist action, and discussed its relation with the EPR spin foam model. In light of the above results, 
we have then re-examined the arguments against the Barrett-Crane model(s), concluding that it can still be 
considered a plausible quantization of 4d gravity as far as the encoding of simplicial geometry is concerned, 
and that further work is needed to either confirm or refute its validity. 
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